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Summary 

Stitched multiaxial warp knit materials have been suggested as viable alternatives to 
laminated prepreg materials for large aircraft structures such as wing skins. Analyses have been 
developed to quantity the effectiveness of stitching for reducing strain energy release rates in 
skin-stiffener debond, lap joint and sandwich debond configurations. Strain energy release 
rates were computed using the virtual crack closure technique. In all configurations, the 
stitches were shown to significantly reduce the strain energy release rate. 

Introduction 

Multiaxial warp knit, 2-D braided, 3-D braided and stitched textile material forms have 
great potential for application to composite aircraft structures and are shown in Figure 1. The 
advantages and limitations of each of the material forms are presented in [1], Among the textile 
material forms that have been developed for aircraft structures, stitched multiaxial waip knit 
materials have proven to be the most viable for large aircraft structures such as wing skins. 

Stitching is expected to suppress the initiation and growth of debonds between the layers 
of the skin and between the stiffeners and skin. The structural skins are typically made of 
between two and ten stacks (layers of material) of 0.055 in. thick carbon warp-knit fabric that is 
layered and stitched with Kevlar yams. The stiffeners and intercostals are fabricated with a 
similar number of stacks of stitched fabric and are then stitched to the skin. Once the preform is 
assembled, resin film infusion (RFI) is used to impregnate the entire structure with epoxy resin 
[1]. Skin-stiffener interface stresses may be large enough to cause a separation between the skin 
and stiffening elements resulting in a delamination or debond. The effects of stitching on 
delamination or debond growth in composites have been examined in [2], This paper 
summarizes the results of fracture mechanics analyses of skin-stiffener debond, lap joint and 
sandwich debond configurations shown in Figure 2. Continuum (2-D and 3-D) and plate finite 
elements are used to model the configuration and the strain energy release rates are calculated 
using the virtual crack closure technique (VCCT) to study the effectiveness of stitching. The 
stitches are modeled as nonlinear spring elements with axial and shear compliances determined 
by experiment. 


Analysis 

A typical stitched wing skin configuration with stiffeners and intercostals is shown in 
Figure 2(a) and a typical stitched sandwich configuration is shown in Figure 2(b). In these 
configurations, debonds can occur between the stiffener and skin, between the layers of skin, or 
between the skin and core. Considerable insight into the behavior of complex configurations 
with debonds can be gained by examining simple configurations such as the skin-stiffener 
debond, lap joint and sandwich debond configurations shown in Figure 3. 

Material and skin thicknesses that are representative of the hybrid 1M7/3501-6 and 
AS4/350I-6 warp-knit fabric stitched composite upper wing skin are considered [1], The 
material consists of IM7 yams in the axial direction (warp direction where warp is also the 
direction of the load) and AS4 yams in the off-axis directions. Each stack of material is 
assumed to be oriented with its primary axis in the v-direction and having a thickness of 0.055 
in. The overall density of the fabric expressed in areal weight is 0.577 oz./ft. 2 , 1.21 oz./ft. 2 and 


0.651 oz./ft. for the forty-five, zero and ninety degree plies, respectively. The equivalent 
laminate stacking sequence of each stack of material is (45/-45/0/9<)/0/-45/45)„ s In these 
analyses, the laminates are assumed to be homogeneous with axial properties determined 
experimentally and all others estimated using the equivalent stacking sequence and classical 
lamination theory as 

Ei pi 1.7 Msi p 12=2.50 Msi u,2=0.40 

E22-5.14 Msi pn=l.77 Msi UipO.30 

E 3 3=1.79 Msi p2.pO.88 Msi o 2 p0.30 

where E\ i, py, u,j (i,j=l,2,3) are the Youngs moduli, shear moduli, and Poissons ratio, 
respectively, and the subscripts 1,2,3 represent the axial (warp), transverse and out-of-plane 
directions, respectively. 

The skin-stiffener debond and lap joint configurations shown in Figures 3(a) and 3(b) were 
modeled using plate elements (STAGS 480, 9-noded quadratic shear deformable, pi ate/shel 1 
elements). The sandwich debond configuration shown in Figure 3(c) was modeled using plane 
strain elements (ABAQUS CPE8, 8-noded plane strain (2-D) quadratic element). Note that the 
offset technique used in plate element modeling [2] is not required for the analyses based on 2- 
D finite elements. For all three debond configurations, the virtual crack closure technique 
(VCCT) is used to calculate strain energy release rates, G, using the nodal forces and 
displacements near the debond front and the debond area created [2], 

Compliance curves for both axial and shear behavior of the stitches were developed in 
reference 3 using flatwise tension and double lap shear tests, respectively. The nonlinear axial 
and shear stiffness, 1 Y axial and K shear , of the spring elements modeling the stitches are extracted 
from these tests. For the purposes of these analyses, the axial and shear responses of the 
stitches are assumed to be independent. A piecewise linear representation of this data is used in 
the finite element model with failure occurring at 58 lb. in tension and 38 lb. in shear [2], 

Results and Discussion 

All three configurations shown in Figure 3 exhibit both mode I and mode II deformations 
at the debond front. I11 the results that follow, both Gi and Gi are shown for the skin-stiffener 
debond and lap joint configurations. Plowever, because of the bimaterial interface between the 
skin and core, only GtoiuI is shown for the sandwich debond configuration. 

The strain energy release rates for the skin-stiffener debond configurations with repeating 
unit boundary conditions (v=0, d,=0 on y=-s y /2) were determined over a range of debond 
lengths in Figure 3(a) with skin and flange thickness, /, of 0.220 inches, and affixed applied 
load of 100 lb./in. In Figure 4, both the mode I and mode II strain energy release rates, G\ and 
Gi, for the unstitched configurations increase with increasing debond length over a range of 
debond lengths, 0ta/ft9.09. The increasing strain energy release rates suggest that a condition 
for unstable debond growth exists in the unstitched configurations. 

Stitching has a significant effect on both Gi and Gi. The mode-I G-value initially increases 
with debond length and then begins to decrease after the debond advances past the first stitch 
(first vertical line in Figure 4), decreasing as additional stitches begin to carry load, eventually 
reaching values of zero for long debonds (u/t>3.41). These zero values for long debonds 
correspond to a region of contact immediately behind the debond front that increases with 
increasing debond length. In contrast, Gi is a weaker function of debond length since the 
shearing displacements are not reduced to zero as the stitches close the debond faces Rather 
the shearing force and displacement are gradually reduced as additional stitches transfer an 


increasing fraction of the load. Thus, even though mode-I is driven to zero for a/t> 3.41 there 
may be sufficient mode-II present to drive the debond for much longer lengths. 

Average debond length vs. load curves were developed from x-ray radiographs of 
unstitched and stitched 27.0 inch long (Z.,) lap joint configurations with a 3.69 inch lap length 
(L 2 ). The thickness, t, of 0.1 10 inches, width, b, of 1.00 inches and stitch spacing, ,v x and ,s„ of 
0.125 inches arc used as shown for the stitched configuration in Figure 3(b). In the 
experiments, debonds grew from only one edge of the unstitched configuration but grew from 
both edges of the stitched configuration. The strain energy release rates w'ere then determined 
using the experimental results of debond length vs. load as input parameters in plate element- 
based finite element analyses. 

Since the single lap joint configuration has a finite width, it may exhibit mode 1 mode II 
and mode III strain energy release rates. Also, although a straight debond front is employed to 
simplify the analyses, a variation in strain energy release rate and stitch force may exist along 
the length (x-) and across the width (>>-) of the lap joint (see Fig. 3(b)). In the results that follow, 
quantities are given at a y-direction location corresponding to: an interior row of stitches near 
the centerline of the stitched configuration (stitch row 1); an exterior row of stitches near the 
edge of the stitched configuration (stitch row 4); and intermediate locations between stitch rows 
1 and 2 and stitch row's 3 and 4. These later locations arc considered only for the stitched 
configuration. They locations denoted as locations A, 3, A and B are located at >=0.0625 in., 
>-0.4j75 in., >=0.125 in. and >=0.375 in., respectively, from the configuration centerline (see 
insert in Figure 5(b)). G-values tor both the unstitched and stitched infinite width 
configurations are also provided. When repeating unit boundary conditions (v=0. 9,=0 on y=0, 
b,2 in Figure 2b) are applied to the models, the calculated G-values are constant in the v- 
direction for the unstitched configuration and periodic in the y-direction for the stitched 
configuration. 


Figure 5(a) shows G-values plotted against debond length for the unstitched configuration 
with a debond at only one end of the single lap joint corresponding to experimentally 
determined debond lengths. Figure 5(a) shows two nonzero components of G (G/ and Gu) near 
the centerline of the finite width lap joint (Location A) and three nonzero components of G (G/, 
G// and Gu/) near the free edge of the finite width lap joint (Location B). Mode I increases at both 
locations as the debond grows. Mode II is the largest component of energy release rate and also 
increases with increasing debond length. Mode 111 has a value near zero (thus, not shown) in the 
interior at location A but increases near the free edge (Location B). 


The energy release rates computed at locations A and B in the finite width configuration 
bound the energy release rates computed for the infinite width configuration (shown as G, and 
Gu in Figure 5(a)). This is a result of the repeating unit boundary conditions in the infinite 
width models preventing anticlastic deformation while the stress free boundary conditions in 
the finite width models do not. Additionally, the repeating unit boundary conditions in the 
infinite width models require mode III to be identically zero along the planes of symmetry. 


Figure 5(b) shows the G-values plotted as functions of debond length for the stitched 
configuration with a debond at each end of the single lap joint corresponding to experimentally 
determined debond lengths. The computed values of Gin Figure 5(b) are at discrete locations 
between the stitch columns and a smooth curve is drawn between the points. Figure 5(b) shows 
one nonzero component of G (G„) at locations A and A and two nonzero components of G (G„ 
and Gw) at locations B and B . For all debond lengths considered at all locations, the mode I 
component is nearly zero. As with the unstitched configuration, mode 11 is the dominant 
component and is : seen to increase with debond length. Although G„ increases with increasing 
debond length, the stitches prevent the debond growth from becoming unstable. As in the 


unstitched configuration, mode III has a value near zero in the interior (Locations A and 4 ) but 
increases near the edge (Locations B and B ). 

The strain energy release rates for sandwich debond configurations modeled with 2-D plane 
strain elements were evaluated for debond lengths, a, between 0.250 in. (a/)y=2.27) and 2.00 in 
(a/'r 18-2) With an applied load, P, of 100 lb. Three core moduli (E con = 5220 psi, E cnn . = 13340 
psi, and E c „n=E /;,<•„ w=l 1.7 Msi) were considered in this study. As shown in Figure 6(a), the 
unstitched configurations exhibit increasing Gr„, a i with increasing debond length, a , for all 
three core moduli considered. Also, for all debond lengths larger values of core modulus 
resulted in decreasing values of G To iai as shown in Figure 6(a). 


The effect of stitching is shown in Figure 6(b) for configurations with four stitches per 
inch. The slope of the G Tota i curve decreases for stitched configurations with foam core and for 
all debond lengths. The decrease is because of the load transfer of the stitches as shown in 
higure 6(b). The slope of the Grmai curve does not begin to decrease for the configuration with 
the stiffest core (F con .—E factsiwer— 1 1.7 Msi) until the third stitch begins to carry load. The G Tot ai 
continues to decrease with both the number of loaded stitches and for larger values of debond 
lengths. The maximum Gw (for the models with the stiffest core) shown in Figure 6(b) is 
0.384 in.-lb./in. 2 at a/tf=4.54 compared with 18.1 in.-lb./in. 2 for the corresponding unstitched 
configuration at a/fy— 18.2. For long debonds (a/tj> 9.09), the strain energy release rate 
approaches near-zero values in the stitched configurations. 


Concluding Remarks 

The effects of stitching on strain energy release rates for skin-stiffener, lap joint and 
sandwich debond configurations were studied using an analysis based on finite elements and 
the virtual crack closure technique. The skin-stiffener debond and lap joint configurations were 
modeled using 9-noded quadratic shear deformable plate/shell elements, while the sandwich 
debond configuration was modeled using 8-noded plane strain quadratic elements. The stitches 
were modeled as discrete nonlinear spring elements with their compliance determined by 
experiment. 


In all configurations studied, the stitches are very effective in the prevention of debond 
growth. The computed values of strain energy release rate (G-values) for all the debond 
configurations studied decrease for longer debond lengths. For the skin-stiffener debond and 
lap joint configurations, the ntode-l G-value (Gi) decreases to near zero values when more than 
three stitches start carrying the load across the debond faces. In contrast, the stitches had less 
effect on mode II and Gn remained significant throughout the range of debond lengths 
considered. Because of the bimaterial interface between the skin and core, only Gw is 
computed for the sandwich debond configuration and it also decreases with increasing debond 
length. 
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